Spatially uniform calibration of a liquid xenon detector at low energies using 83m Kr 
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A difficult task with many particle detectors focusing on interactions below ~100 keV is to perform 
a calibration in the appropriate energy range that adequately probes all regions of the detector. 
Because detector response can vary greatly in various locations within the device, a spatially uniform 
calibration is important. We present a new method for calibration of liquid xenon (LXe) detectors, 
using the short-lived 83m Kr. This source has transitions at 9.4 and 32.1 keV, and as a noble gas like 
Xe, it disperses uniformly in all regions of the detector. Even for low source activities, the existence 
of the two transitions provides a method of identifying the decays that is free of background. We 
find that at decreasing energies, the LXe light yield increases, while the amount of electric field 
quenching is diminished. Additionally, we show that if any long-lived radioactive backgrounds are 
introduced by this method, they will present less than 67 x 10 -6 events/kg/day in the next generation 
of LXe dark matter direct detection searches. 

PACS numbers: 29.40.Mc; 78.70.-g; 61.25.Bi 



I. INTRODUCTION 

Liquid xenon (LXe) particle detectors are commonly 
used in the fields of dark matter direct detection [1-5] , 
neutrinoless double-beta decay searches [6] , collider ex- 
periments [7] , and in nuclear medicine [8, 9]. In these 
applications, energy deposition is measured by count- 
ing scintillation photons, ionization, or both. The high 
electron mobility [10] and scintillation properties [11] of 
this liquid contribute to its versatility. LXe response de- 
pends on the electronic stopping power (dE/dx) to the 
incoming particle, which itself depends on both the iden- 
tity and energy of the particle. It is therefore neces- 
sary to calibrate a detector's energy scale with a source 
whose response is known relative to the particles under 
study. One common such "reference source" is 57 Co, 
which emits y-rays predominantly at 122 keV. 

Dark matter direct detection experiments search fol- 
low energy nuclear recoils caused by the scattering 
of weakly interacting massive particles (WIMPs) with 
atomic nuclei [12-14]. There are two main problems 
involved in using 57 Co to calibrate LXe detectors for 
this application. The first is that the y-ray energy is 
much higher than the recoiling nuclei energy produced 
by WIMP interactions. Second, the attenuation length 
of 122 keV y-rays in LXe is ~2.5mm, and hence the en- 
ergy deposition will be highly localized as compared with 
the tens of cm typically characterizing the size of such de- 
tectors. The two problems are actually compounded, be- 
cause the attenuation length of y-rays decreases as their 
energy decreases, and therefore sources providing lower- 
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energy y-rays will give an even more localized response 
than 57 Co. The topic of localization is an issue for point 
sources placed outside the detector, but also for point 
sources placed inside the detector. In the latter case, the 
source must be attached to a mounting device of some 
kind; for low energy y-ray sources, any device used for 
this purpose will likely block some of the scintillation 
light and potentially distort any existing applied electric 
fields. It is therefore not possible to calibrate a detec- 
tor with an internal point source under the same condi- 
tions that would exist during the actual measurement. To 
avoid these difficulties, short-lived noble gas sources can 
be used which diffuse uniformly in LXe. The XENON10 
experiment used the metastable 131m Xe [4, 15]. This 
source solves the second problem (spatial localization), 
but its 164 keV transition does not overcome the prob- 
lem of an appropriate energy scale. Additionally, due to 
its half-life of twelve days, the detector must sit for ap- 
proximately 2.5 months following a calibration until the 
source activity has dropped to 1% of its initial value. 

A promising alternative solution is to use the 
metastable 83m Kr, first proposed in [16]. This source 
has been used for calibrations of detectors in the Large 
Electron-Positron Collider [17, 18], as well as in the KA- 
TRIN experiment which attempts to measure the elec- 
tron neutrino absolute mass [19]. 83m Kr should diffuse 
uniformly in a LXe detector, addressing the issue of 
spatial uniformity (while perfect uniformity of deposi- 
tion is not explicitly required, a uniform deposition al- 
lows for calibration of the full detector in the minimum 
time). Additionally, its two de-excitation lines at 9.4 and 
32.1 keV lie in the energy range of interest for dark matter 
direct searches, and its half-life of only 1.8 hours allows 
for a short turnaround time following measurement. Si- 
multaneous with the work presented here, another group 
has performed a demonstration of a similar 83m Kr in- 
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troduction method in a single-phase LXe chamber [20] . 
Here we present a successful implementation of this cal- 
ibration source in a small dual-phase LXe prototype de- 
tector. Furthermore, we show results of measurements of 
the LXe energy scale linearity, evolution of energy reso- 
lution with energy, effects of LXe response under applied 
electric fields, and set limits on the level of long-lived 
radiocontaminants introduced by this method. 



details on the production process, we refer to [22]. In ad- 
dition to 83 Rb, 84 Rb and 86 Rb isotopes are also produced, 
however, they decay to stable Kr isotopes and hence in- 
troduce no radioactive backgrounds. Since 83 Rb decays 
with a half-life of 86.2 days, the source strength decreased 
to ~3kBq by the end of these measurements. 



B. Detector System 



II. EXPERIMENTAL METHODS 

A. The 83m Kr source 



83m Kr is produced by the decay of 83 Rb via pure elec- 
tron capture. Following this process, the nucleus is left in 
any of 83 Kr's many excited states lying below the Q- value 
of the Rb decay (910 keV). Regardless of the initial kryp- 
ton excited state, the nucleus rapidly de-excites within 
picoseconds to the isomeric state 83m Kr, located 41.5 keV 
above the ground state. Isomeric krypton decays with a 
half-life of 1.83 h to the first 83 Kr excited state (9.4 keV), 
which then decays to the ground state with a half-life 
of 154 ns [21]. The decay scheme of 83m Kr is shown in 
Figure 1 , indicating that most of the released energy is 
carried by internal conversion and Auger electrons [18]. 
The 6kBq 83 Rb source used in this study was produced 
at the Nuclear Physics Institute, Rez (Czech Republic). 
This institute also provides 83 Rb for the KATRIN ex- 
periment [19]. The parent 83 Rb is produced in the U- 
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The 'Xiirich' detector, used for studying the 83m Kr de- 
cay is shown schematically in Figure 2. It is a dual-phase 
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FIG. 2: Schematic of the dual-phase Xiirich detector. The 
PTFE structure holds the PMTs and grid electrodes (see 
text), defining an active region 3.5cm in diameter and 3cm 
high. 



FIG. 1: The decay scheme and branching ratios of Kr. The 
decay always passes through two transitions, giving mostly 
internal conversion (IC) and Auger (A) electrons. A small 
amount of the energy is carried by gamma- (7) and X-rays 
(X) [18]. 

120M cyclotron from the reaction nat Kr(p,xn) 83 Rb by ir- 
radiating a medium-pressure krypton target with 27McV 
protons. The product, deposited on the target chamber 
walls, is then washed into several tens of milliliters of high 
purity water (<0.07/iS/cm). An appropriate amount of 
the target eluate is then absorbed in zeolite beads (2 mm 
diameter, Merck), which acts as a molecular sieve. Ze- 
olite was chosen due to its ability to allow for efficient 
emanation of 83m Kr in vacuum, while exhibiting high re- 
tention of the mother 83 Rb in its porous structure. For 



(liquid and gas) LXe time projection chamber (TPC). 
The system was built at the Universitat Zurich specif- 
ically for testing liquid xenon's response to low energy 
interactions. The stainless steel vessel in Figure 2 is lo- 
cated inside a vacuum cryostat, with cooling provided 
via a copper cold finger immersed in liquid nitrogen. 
The temperature and pressure are held constant at 175 K 
and 1.8 bar, respectively, and the detector was operated 
stably for approximately 2.5 months. Xiirich's cylindri- 
cal active region, 3.5 cm in diameter and 3 cm in height 
(80.8 g of LXe), is defined by a polytetrafluorocthylcne 
(PTFE) cylinder on the perimeter and grid electrodes 
above (gate) and below (cathode). A third grid electrode 
(anode) is located above the gate grid, with the liquid 
level lying between the gate and anode grids. Two Hama- 
matsu R9869 [23] photomultiplier tubes (PMTs) view the 
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active volume, one from below and one from above. A 
total of 1.76 kg is used to fill the stainless steel vessel. 

The cathode and gate grids apply an electric field of 
typically ~ 1 kV/cm which is used to drift electrons away 
from an interaction site towards the gate grid. Once the 
electrons pass through the gate grid, they arrive at the 
liquid surface and are extracted to the gas by an electric 
field of ~ lOkV/cm that then accelerates the electrons 
through the gas until they collect on the anode grid. The 
high voltage applied to the grids is supplied by a CAEN 
model A1526 module. During their transit through the 
gas, the electrons collide with Xe atoms with sufficient 
energy to produce scintillation light. Therefore, the typ- 
ical result of a particle interaction is a prompt scintilla- 
tion signal (SI) emitted from the interaction site itself, 
followed by a delayed scintillation signal (S2) produced as 
the electrons travel through the gas under the influence 
of the extraction field. In this way, both the scintilla- 
tion and ionization signals are measured by the PMTs. 
This technique is used for charge readout because it pro- 
vides superior amplification over more traditional meth- 
ods [24, 25]. Additionally, the z-position of the event can 
be inferred by the delay time between the SI and S2 sig- 
nals since the electron drift velocity is well known as a 
function of applied field [10]. 

The light yield is defined as the number of photo- 
electrons (p.e.) emitted from the PMT photocathodes 
per unit energy, and is customarily quoted based on 
the primary emission of 57 Co. When Xiirich is oper- 
ated in single-phase mode with the liquid level above 
the top PMT, the 57 Co source produces ^10p.e./keV. 
In the dual-phase mode used in this study, where the 
liquid level lies below the top PMT, some of the scintil- 
lation light is refracted and reflected at the liquid sur- 
face. Though some of the reflected photons may be de- 
tected by the bottom PMT, roughly 35% are lost. The 
result is a significantly larger SI signal in the bottom 
PMT compared to the top (70% on bottom, 30% on 
top), and an overall reduced light yield as compared 
with the value taken in single-phase mode. We mea- 
sure the dual-phase 57 Co zero-field light yield to be 
6.38±0.05(stat)±0.36(sys)p.e./keV, with 11.5% resolu- 
tion (a/n). The systematic uncertainty is taken from the 
fluctuations of this light yield throughout the run, and 
the statistical uncertainty is the combination of fit uncer- 
tainties from each 57 Co zero-field data set. The spectrum 
obtained from one 57 Co calibration is shown in Figure 3. 

The PMTs were produced with photocathodes specif- 
ically designed to have their peak quantum efficiency 
(^25% bottom, ^35% top) at 178 nm, the wavelength 
of Xe scintillation light [26] . Negative high voltage of ap- 
proximately -900 V, supplied by a NHQ 225M high volt- 
age module, is distributed to the photocathodes and 12 
dynodes by voltage dividers made in-house on a PTFE 
substrate with a total impedance of ^140 Mf2. Each 
PMT is calibrated in situ by a pulsed blue LED, period- 
ically and before each measurement, to obtain the single 
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FIG. 3: Spectrum obtained from 57 Co at zero field and Gaus- 
sian fit. The run-averaged light yield is 6.38 p.e. /keV. 



p.e. response, and hence the PMT gain (~4xl0 6 at these 
voltages). The PMT signals are fed to an external fast 
voltage amplifier (Phillips 777) having two outputs; one 
output is connected directly to the ADC (Acqiris model 
DC436 lOOMS/s), while the other is fed to the trigger- 
ing system. The triggering system is set to discriminate 
each PMT individually at -~1.5 p.e., requiring coincidence 
between the two PMT channels. The efficiency of this 
trigger, determined by Monte Carlo method, is 50% at 
1.6 keV, rising to 95% at 4keV and 99% at 5.4 keV, well 
below the 32.1 keV events that provide the trigger in a 
83m Kr decay. In the offline signal processing, the pulse 
areas are converted to a value in number of p.e. by using 
the measured PMT gains. 

In order to maintain the purity of LXe at a level suffi- 
cient for optimal light and charge collection, the xenon is 
constantly purified with a hot metal getter, SAES Mono 
Torr model PS3MT3R1. The purification is handled by 
the gas system shown schematically in Figure 4. Liquid 
is drawn up from the detector by a closed diaphragm re- 
circulation pump, Enomoto MX-808ST-S, where it is va- 
porized upon leaving the cryostat. After passing through 
the pump, it passes through the getter and back into the 
Xiirich detector where it is recondensed. In this way, the 
xenon is constantly purified in a closed recirculation loop. 

83m Kr is introduced into the flow of the closed recircu- 
lation circuit by means of a single port with a valve. The 
zeolite beads containing the 83 Rb reside in a small cham- 
ber filled with the same xenon gas that flows in the gas 
system. Gaseous 83m Kr emanating from the 83 Rb de- 
cay may then diffuse into the recirculation circuit, its 
introduction being easily controlled by either opening or 
closing the valve at the port, denoted as the Rb valve. 
Due to the rather long half-life of 83 Rb (86.2 d), it is im- 
perative that no trace of this mother radionuclide enters 
the system if it is to be used in a low background ex- 
periment. Rb might potentially enter the system by one 
of two ways: as a vapor, which is very unlikely since its 
volatility under common laboratory temperatures even 
under vacuum is not significant; or as an aerosol formed 
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FIG. 4: Schematic of the gas system. The recirculation pump 
draws liquid up through the recirculation siphon tube in the 
LXe chamber where it is vaporized while leaving the cryostat 
space, passed through a getter for purification, and the re- 
condensed inside the chamber. The zeolite beads containing 
S3 Rb are located inside the chamber labeled 'Rb' and expo- 
sure to the gas system is controlled by the valve labeled 'Rb 
valve'. 



from small particles of the zeolite itself. Aerosol break- 
through is not entirely excluded, and therefore a 0.5 /Ltm 
aerosol filter is placed between the Rb chamber and the 
Rb valve in order to prevent any 83 Rb from entering the 
recirculation loop. Measurements done to assess the level 
of 83 Rb introduced in the the system are discussed in sec- 
tions III and IV. 



III. ANALYSIS AND RESULTS 

Once the 83m Kr has entered the LXe, a 32.1 keV tran- 
sition might occur in the active region, which will then 
be followed by the 9.4 keV transition. A 83m Kr decay 
is, therefore, indicated by two SI pulses whose separa- 
tion in time is characterized by a decaying exponential 
with i 1 / 2 =154ns. Some of these transitions will occur 
too close in time to be resolved separately, giving a sin- 
gle 41.5 keV pulse; however, the strength of this signal is 
well below the background level in the Xiirich detector. 
On the other hand, many of the 83m Kr decays have a 
double SI structure, while only a small fraction of non- 
83m Kr decay events share this feature. An example of 
the PMT response from a 83m Kr decay is seen in Figure 
5 (top). 

The events with such a double SI structure are shown 
from one data set in Figure 5 (bottom), with the area of 
the first pulse plotted versus the area of the second pulse. 
In this space, it is evident that the 83m Kr decays form a 
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FIG. 5: (Top) PMT output from a 83m Kr decay. In this dou- 
ble pulse of primary scintillation light (SI), the first pulse cor- 
responds to the 32.1 keV transition with the second pulse re- 
sulting from the 9.4 keV transition. (Bottom) The area of the 
first SI pulse versus the area of the second, for events show- 
ing this characteristic two-pulse structure. Shown are distri- 
butions taken before Rb exposure ('Background') and during 
Rb exposure (' 83m Kr'), demonstrating that the population of 
83m Kr decays is clearly separated from background events. 
The box represents the energy cuts used as the 83m Kr accep- 
tance window. 



population of events that is clearly separated from back- 
ground. The box indicates the energy cuts for first and 
second SI pulses used to identify 83m Kr decays; before 
opening the Rb valve, background data show no events 
within this box. After the Rb valve has been opened, the 
rate of 83m Kr decays in the total LXe volume increases to 
the 20 Bq level in roughly 10 h. In order to further check 
that these are indeed 83m Kr decays, the distribution of 
SI delay times (i.e. the time between the first and sec- 
ond SI pulses), Atgi, of events within the box of Figure 
5 (bottom) is fit with a decaying exponential. The result 
of the fit, shown in Figure 6 (top), gives ti/ 2 = 156±5ns, 
consistent with the published value of 154.4±l.lns [21]. 
This excellent agreement validates the claim that these 
events are indeed caused by 83m Kr decays. 

Due to the shaping of the PMT signals by the various 
DAQ components, multiple SI pulses that are delayed 
by less than ^100 ns cannot be separately resolved. Ad- 
ditionally, the signal is required to be 'clean' (i.e. flat 
baseline) two samples before and after the pulse, in 
order to register as a positive SI identification during 
the offline processing of the data. This makes the effi- 
ciency for detecting multiple SI pulses less than unity 
for Aisi < 250ns, as is obvious from Figure 6 (top). 
Therefore, the double SI selection cut detects 83m Kr de- 
cays with an efficiency of approximately 32% under these 
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FIG. 6: (Top) The distribution of delay times between first 
and second SI pulses for events in the 83m Kr acceptance win- 
dow. An exponential fit to the distribution gives a half-life 
of 156 ± 5 ns, consistent with the published value of 154 ns. 
(Bottom, color online) Spectra from the two 83m Kr transi- 
tions, summed over all runs taken at zero field. 



conditions. 

The spectra, in p.e., obtained at zero field from the 
two transitions of 83m Kr are displayed in Figure 6 (bot- 
tom). A Gaussian function is fit to each spectrum that is 
used to determine the light yield and energy resolution, 
shown in Tabic I. As mentioned in section I, 57 Co emits 
primarily 122 keV 7-rays. However, there is a small ad- 
ditional contribution from 136 kcV. The two lines, how- 
ever, are not distinguishable from one another due to 
the detector's energy resolution and instead give a single 
peak, whose weighted average energy is 123.6 kcV. The 
measurements suggest a rise in the light yield at lower 
energies, consistent with behavior previously observed in 
LXe[27] and also in the XENON10 detector [28] . The 
peak resolutions (a / ' (i) are also shown at zero field. 

As mentioned in subsection II B, most LXe detectors 
use an applied external electric field in order to collect 
electrons emitted from the interaction site. As the ap- 
plied field is increased, more and more electrons leave 
the interaction, suppressing the recombination process 
that contributes photons to the scintillation signal. The 
result is that both the scintillation and ionization re- 
sponses vary strongly with applied field, with the two 
signals exhibiting anti-correlation. It is then crucial that 
the scintillation dependence on the applied field, called 
field quenching, be known quantitatively for any calibra- 
tion sources. Figure 7 shows the light yield as a function 
of the applied field, normalized to the zero field value, of 
the two 83m Kr transitions and the 57 Co line. The uncer- 



tainty in the light yield is dominated by a 2% systematic 
uncertainty taken from the measured fluctuations in the 
PMT gain over the duration of the run. The horizon- 
tal positions are determined by electrostatic field simula- 
tions of the detector in each voltage configuration used; 
horizontal uncertainties are the 1-cr variation of the field 
over the active volume. The simulations were carried 
out using the COMSOL simulation package (commer- 
cially available) [29], and verified with software written 
in-house. 

The time scale of the ionization signal, 1-2 (is, does 
not permit the two 83m Kr transitions to be resolved sep- 
arately, and instead the S2 signal contains the combina- 
tion of charge emitted from both decays. This 41.5 keV 
summed-signal ionization yield is also shown in Figure 
7 normalized to Qq, the theoretical total amount of ini- 
tial charge produced prior to electron-ion recombination. 
This value is determined by plotting the SI peak posi- 
tions versus the S2 peak positions from data taken at var- 
ious applied fields, shown in Figure 8. As SI and S2 are 
anti-correlatcd, these data lie along a line having nega- 
tive slope, with the line's intercepts representing the total 
number of quanta, ions plus excitons (Ni on + N ex ). For 
electronic recoils, the ratio of excitons to ions, N ex /Ni on , 
is taken to be 0.06 [30], and hence Qo is 94.3% the value 
of the S2 intercept. The horizontal positions and error 
bars are determined in the same manner as those of the 
scintillation yield measurements, while the vertical error 
bars are instead dominated by the statistical errors in the 
peak fits and the uncertainty in Qq. 
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FIG. 7: (Color online) Field quenching, defined as the light 
yield of a spectral line divided by the light yield obtained 
at zero field, or S(E)/S(0). Data collected from 57 Co (open 
black squares) are consistent with those reported in [31] (solid 
grey diamonds). Dashed lines correspond to a fit parameteri- 
zation described in the text. Also shown is the field-dependent 
charge collection of the combination of both 83m Kr transi- 
tions, Q(E)/Qo; the two transitions occur too close in time 
for their ionization signals to be individually resolved. 
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TABLE I: The measured zero-field light yield (L.Y.), peak resolution (Res.), and field dependence fit parameters, a;. The row 
following 41.5 keV gives the charge collection of the summed signal. Uncertainties shown in light yield are statistical only; 
because these two peaks are taken from identical events, their systematic uncertainties are highly correlated, and hence do not 
affect the significance of the relative rise in light yield. 



E (keV) 


L.Y.(p.e./keV) 


Res. (<t//u) 


ai 


a 2 (lO^cm/V) 


a-i 


9.4 


6.74±0.06 


20.0% 


-0.35±0.06 


6.3±3.0 


1 


32.1 


6.43±0.04 


14.4% 


-0.55±0.03 


8.9±1.6 


1 


41.5 






0.406±0.006 


17±2 


0.074±0.012 


123.6 


6.38±0.05 


11.5% 


-0.679±0.007 


12.6±0.5 


1 



The data are fit with a three-parameter function based 
on the Thomas-Imcl box model for electron-ion recombi- 
nation [32] , given by, 



S(E) Q(E) 
S(0)' Q 



= CL1CL2E In 1 



1 



a 2 E 



03, 



(1) 



where E is the electric field strength, and S, Q are 
the scintillation and ionization yields, respectively. This 
model is used only to provide a convenient parameter- 
ization of the data, and not to infer fundamental LXe 
physical properties from the results of the fits. The ai 
are the parameters of the fit, shown in Table I. Because 
the scintillation yields are normalized to the value at zero 
field, CI3 is unity and the function therefore contains only 
two free parameters for the scintillation yield data. At 
decreasing energies, we observe a consistent decrease in 
the level of field quenching. 
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FIG. 8: (Color online) The peak position in S2 versus SI space 
for the 41.5 keV emission of 83m Kr. The data are taken from 
applied fields ranging from lOOV/cm to lkV/cm. The line 
is fit to the data having vertical and horizontal intercepts Is2 
and 7gi, respectively; these intercepts indicate the location of 
N ion +N ex . 



The energy of an event can also be measured by count- 
ing the total number of quanta, Ni on -\~N ex . This is called 



the combined energy scale (CES), and is constructed by 
forming a linear combination of the scintillation and ion- 
ization signals, aSl + /3S2, such that n 7 = aSl and 
n e = /3S2, where n 7 and n e are the number of emit- 
ted photons and electrons, respectively. The coefficients 
a and (3 can be found from the plot of SI versus S2 (Fig- 
ure 8), by a = E/{WI S1 ) and fi = E/(WI S2 ), where 
E is the deposited energy, 11^=13.5 eV is the average en- 
ergy required to produce a single quanta (electron or pho- 
ton) [33], and isi(S2) is the S1(S2) intercept in units of 
p.e.. The CES has the advantage that it is not affected 
by correlated recombination fluctuations which dominate 
the SI resolution over most energies [32], and hence gives 
an energy estimate with better resolution than SI or S2 
alone. For example, the Sl-only and S2-only peak resolu- 
tions of the 41.5 keV peak taken at 500V/cm are 14.2% 
and 20.1%, respectively, while the resolution of the CES 
peak at this field is 10.0%. 

The delay time between SI and S2 gives the drift time 
of the electrons, and hence the z-position of the inter- 
action. One important motivation for using this source 
is that it should disperse uniformly in the active LXe 
volume, providing a spatially-uniform calibration. The 
summed z-position distribution of 83m Kr events taken 
at drift fields from 100-1000V/cm is shown in Figure 9 
(top). The observed z-dependent rate is flat with varia- 
tions consistent with statistical fluctuations on each bin. 
With this uniform calibration, the position-dependence 
of the detector's response can be measured and cor- 
rected for. Most of the SI signal is detected by the bot- 
tom PMT, and therefore one expects to see a light yield 
that is a monotonically decreasing function of z-position 
(i.e. more light is collected from events occurring close to 
the bottom PMT than for events close to the top). Fig- 
ure 9 (bottom) shows the light yield of the 83m Kr decays 
at all positions along the z-axis between the cathode and 
gate grids for the data run at 1 kV/cm; solid lines are the 
band ccntroids, shaded bands cover ±1<t. In both tran- 
sitions, the light yield at the cathode (bottom of active 
region) is a factor of 1.3 higher than the light yield at the 
gate grid (top of active region). Data collected at lower 
applied fields shows consistency with this behavior. 

Although the 83m Kr decays away in a matter of hours, 
the 83 Rb will live for nearly 1.5 yr before decaying below 
1% of the initial activity. If this technique is to be used in 
low-background experiments, it is then imperative that 
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FIG. 9: (Top) Rate of 83m Kr decays as a function of 
z-position, indicating a uniform concentration. (Bottom, 
color online) Measured z-dependence on the light yield from 
83m Kr's two transitions taken at lkV/cm. The solid lines 
indicated the band centers, with ±lcr covered by the shaded 
areas. Both lines show a light yield at the cathode that is a 
factor of 1.3 larger than at the gate grid. 



no 83 Rb atoms enter the system, and instead must remain 
trapped within the zeolite or the filter. In order to test 
this, the valve to the 83 Rb chamber was closed. The rate 
of 83m Kr decays is expected to decrease exponentially 
to zero during the following day; however, if 83 Rb has 
entered the system, the rate vs. time will behave as an 
exponential decay with a vertical offset. No such offset 
was observed in the 83m Kr rate following the closing of 
the Rb valve. Indeed, 2.5 h of data collected one day after 
closing the Rb valve resulted in zero observed events. 
This null observation corresponds to a one-sided Poisson 
upper limit of <800/zBq (90% C.L.) of residual 83m Kr in 
the active region. 



IV. DISCUSSION 

Testing the low energy response of LXe is generally 
rather difficult, and therefore 83m Kr provides a unique 
tool for such measurements. The rate of 83m Kr decays 
studied in this work is quite low as compared with the 
rate of background events due to natural radioactivity 
and cosmic rays, contributing roughly 0.5% of the total 
triggers. However, the double SI structure of these de- 
cays, and energy cuts used, enable their measurement in a 
virtually background-free regime. It is therefore not nec- 



essary to use a low background setup simply to study this 
weak source. Simultaneous with this work, a demonstra- 
tion of 83m Kr introduction to a single-phase LXe cham- 
ber by a similar technique has been performed by another 
group [20]. 

The light yield and energy resolution at low energies 
are of particular relevance for dark matter direct detec- 
tion searches. As indicated in Table I, the light yield 
increases at low energies. Although an accurate quan- 
titative understanding of this process is incomplete, the 
observed behavior can be understood qualitatively in the 
following manner. The electronic stopping power of elec- 
trons in LXe increases at decreasing energies [34] , and 
thus the ionization density produced by a recoiling elec- 
tron increases along the track, with the highest densities 
concentrated at the track's end. Because of this, the 
overall ionization density caused by a low energy recoil- 
ing electron will be greater than for an electron of higher 
energy. The electrons and ions produced along the track 
will rapidly recombine and produce scintillation photons 
as they de-excite to their ground states. The strength of 
recombination is correlated with the ionization density, 
because the characteristic electron-ion distance is shorter 
for higher ionization densities. Even at zero applied elec- 
tric field, not all of the electron-ion pairs produced by 
a recoiling particle will recombine to give scintillation 
photons [II]. It is then expected that the zero field re- 
combination is stronger at lower energies (higher dE/dx), 
giving a higher overall light yield. 

In [22], 4 h baking at 300 °C resulted in no observed loss 
of 83 Rb. The observation in this study of no 83m Kr de- 
cays, one day following the closing of the Rb valve, sets 
an upper limit of 800 ^Bq of residual 83m Kr inside the 
active region of the Xiirich detector. Prior to this, the 
Rb valve had been opened for a total of 150 hours during 
the run. The risk of Rb contamination increases with the 
amount of time that the valve is opened, and so this up- 
per limit can be normalized to exposure time (150 h). 
Moreover, since the source is exposed to the gas sys- 
tem (and not the detector), the total activity in the LXe 
chamber should be independent of the detector size, and 
should instead depend on only flow rate and method of 
deployment. The limit of 800 /iBq in the active region 
(0.08 kg) can be scaled to the total amount of LXe in the 
chamber (1.76kg), and normalized to the exposure, to 
give <120/iBq/h of residual 83m Kr in the whole liquid 
volume (assuming the 83m Kr concentration outside the 
active region is uniform and equal to the concentration 
inside the active region). The branching ratio of 83 Rb to 
83m Kr is 75%, which means this limit on residual 83m Kr 
is a limit of <160/iBq/h of residual 83 Rb contamination. 

To understand how this upper limit would affect an 
actual dark matter search, a 300 kg detector with 100 kg 
fiducial mass is taken as an example. A detector of this 
size is typical of the proposed next generation of LXe 
dark matter searches [35, 36]. An exposure to 83 Rb of 
10 h would be sufficient, under these conditions, to pro- 
vide adequate statistics for such a calibration (~1000 
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Kr events/kg). Our upper limit of 83 Rb contamination 
translates to a residual rate of <0.46decays/kg/day in 
this 300 kg. Even if this amount of 83 Rb was present 
in the system, the vast majority of decays would not 
introduce dangerous backgrounds. In order for a back- 
ground event to be 'dangerous' (i.e. appear in the WIMP 
signal acceptance window), it must have two features: 
(1) it must produce a single scatter event; (2) the event 
must deposit a small amount of energy that is within the 
WIMP search energy window. An additional feature that 
dual-phase LXe TPCs like Xiirich have is the ability to 
reject electronic recoils on an event-by-event basis at the 
level of ~99.9% based on the ratio S2/S1 [37]. However, 
statistical fluctuations can cause a small fraction of elec- 
tronic recoil events to yield a S2/S1 ratio similar to values 
characteristic of a nuclear recoil from WIMPs, and thus 
the overall background level must be minimized as much 
as possible. Any 83m Kr decays in the active volume would 
not present a problem because they would either have a 
double SI structure (and could be vetoed on that basis), 
or would give 41.5 keV, outside of the WIMP search re- 
gion. The only possibility for a dangerous background 
is from one of the 7-rays produced as the initial excited 
83 Kr decays to the metastable state. These 7-rays are 
mostly emitted in the range of 500-600 keV; again, to be 
dangerous they arc required to single-scatter in the fidu- 
cial region, which is highly unlikely given their 3-4 cm 
attenuation length. With 83 Rb contamination at the 
level of our upper limit, Monte Carlo simulations indi- 
cate that 0.46 decays/kg/day would contribute less than 
67 /iDRU of single scatters in the WIMP search energy re- 
gion (1 DRU=1 event kg -1 day -1 keV -1 ). The projected 
7 background rate in [35] and [36] due to natural radioac- 
tivity in the detector materials alone is roughly 1 mDRU, 
fully fifteen times greater than our upper limit on the 
83 Rb background. 

After a calibration with 83m Kr, the stable 83 Kr will 
remain in the system indefinitely unless some action is 
taken to specifically remove it. However, the amount of 
Kr remaining from a 10 h exposure as described above 
will be miniscule; less than 10 6 atoms total, which cor- 
responds to a concentration of roughly 1 part in 10 21 for 
300 kg of Xe. Even if this remaining concentration was 
higher, Kr will not adversely affect detector functions. 
The transport of electrons through the Xe will not be 
diminished since Kr is chemically similar to Xe. Addi- 
tionally, Kr does not absorb Xe scintillation light [38] and 
therefore will not impede light collection. 

In addition to 83 Rb contamination, water and oxygen 
trapped in the zeolite might also enter the system. While 
these elements do not pose a problem in the context of 
radioactive backgrounds, they could affect the charge col- 
lection and light yield. Before the Rb valve was initially 
opened, the Rb chamber was evacuated to the level of 
10~ 6 mbar with a turbomolecular pump at room temper- 
ature. The Rb valve was then open continuously for ap- 
proximately four days, following which diminished charge 
collection was observed. The Rb valve was then closed 



and the purification system allowed time to restore the 
LXe purity to a level adequate for negligible charge loss. 
In subsequent measurements, the Rb valve was toggled in 
cycles of 20 h open, 4 h closed, with charge collection pe- 
riodically monitored; no charge loss was measured under 
these conditions. It is likely that the impurity content in 
the zeolite was depleted in the initial four days of expo- 
sure, and had left the system by the time the cycles of 
20 h exposure began. In a subsequent run, the Rb cham- 
ber was baked for 24 hours at 120 °C prior to exposure, 
following which no effect on the electron lifetime was ob- 
served. At no time was any effect on the light collection 
seen. 

The Xiirich detector will be used in the future to study 
light and charge collection in LXe in response to nuclear 
recoils, scintillation pulse shape studies, and in studies 
related to research and development leading to a ton- 
scale LXe dark matter search. 



V. SUMMARY 

This work presents a method to utilize the radioactive 
source 83m Kr in a LXe detector. The two transitions of 
the 83m Kr decay produce two primary scintillation pulses, 
which allow these events to be clearly separated from 
background. Additionally, the dual-phase detector used 
to study this source shows a spatially uniform distribu- 
tion of decays, validating one of the main motivations for 
studying this source. The data collected with this source 
show an increase in the light yield at decreasing energies, 
indicating a corresponding increasing level of zero-field 
electron-ion recombination. A corresponding weaker field 
dependence of the light yield at decreasing energies is also 
observed. An alternative energy scale based on the com- 
bination of scintillation and ionization signals was shown 
to provide superior energy resolution at 41.5 keV com- 
pared with scintillation or ionization alone. 

The technique used to introduce 83m Kr involves ex- 
posing the system to 83 Rb, the presence of which in a 
low-background experiment could yield serious problems 
if contamination occurs. This measurement provides an 
upper limit on the level of 83 Rb contamination intro- 
duced to the system of <160/iBq/h exposure. In the 
context of the next generation of LXe dark matter direct 
detection searches, this upper limit translates to a max- 
imum introduced background of 67 /iDRU, far below the 
1 mDRU of intrinsic backgrounds. 
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